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We show that, under rather general assumptions, the phase diagram of a quasi-one-dimensional 
repulsive Fermi system consists of two ordered phases: the density wave, spin or charge, and the 
superconductivity. It is demonstrated that the symmetry of the superconducting order parame- 
ter is a non-universal property sensitive to microscopic details of the model. Three potentially 
stable superconducting states are identified: they are triplet /-wave, singlet d^2_y2-wa.ve, and dxy- 
wave. Presence of multiple competing superconducting states implies that for a real material this 
symmetry is difficult to predict theoretically and hard to probe experimentally, since artifacts of 
theoretical approximations or variations in experimental conditions could tip the balance between 
the superconducting phases. 



I. INTRODUCTION 

Many theorists hold that the superconductivity in 
quasi-one-dimensional (QID) materials may be explained 
with the help non-phonon many-body mechanism [1, 2, 
3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17]. Yet, 
there is substantial disagreement about the details of 
such mechanism. Moreover, different theories predict 
different order parameter symmetries: p—, d.j.2_y2-, f- 
, or da^y-wave [8, 9, 10, 11]. The experimental findings 
[18, 19, 20, 21, 22, 23, 24, 25] have been unable to resolve 
these controversies. 

In this article we explain why establishing the symme- 
try of the order parameter is such a hard task. It will 
be demonstrated that within the framework of the non- 
phonon mechanism there are three metastable supercon- 
ducting states: dx2_y2-, /-, and dxy-wetve, which compete 
to become the ground state. This circumstance has an 
important implication for both theory and experiment: 
to find the ground state symmetry one has to resolve 
small energy differences between the contesting states. 
Thus, an approximate theoretical scheme might calcu- 
late these differences completely wrong; a slight variation 
of the model may cause a phase transition. An experi- 
mentalist has to keep in mind that, despite the obvious 
similarities between various families of the QID super- 
conducting materials in chemistry and structure, despite 
"universality" of their phase diagram [25] , the supercon- 
ducting order does not have to be identical in all these 
metals. Moreover, even an individual sample may ex- 
perience phase transitions between different types of su- 
perconductivity as the pressure or magnetic field change 
(e.g., [24]). This suggests that experimental detection 
of the order parameter symmetry should be done for a 
specific material, rather than a class of materials; a re- 
liable theoretical prediction of this symmetry is nearly 
impossible. 

Let us briefly outline main ideas of the discussion. Due 



to presence of strong many-body effects, the QID metal 
cannot be treated by the usual mean-field theory. Thus, 
to circumvent this difficulty, we derive the low-energy 
effective Hamiltonian. Unlike the original microscopic 
model, the anisotropy of the effective model is low, and 
it may be studied with the help of the mean-field ap- 
proximation. Application of mean-field theory to the ef- 
fective model [1, 2] allows us to find the phase diagram 
of the QID metal. At good nesting the system freezes 
into a spin or charge density wave. The antinesting de- 
stroys the density wave phase and makes superconduc- 
tivity possible. In superconducting phase we find three 
order parameters, which may be stable in our model. 

The paper is organized as follows. The QID model 
and its effective low-energy description is presented in 
Sect. II. The phase diagram is obtained in Sect. III. 
Sect. IV is reserved for discussion of the results. 



II. QID MODEL AND ITS EFFECTIVE 
HAMILTONIAN 

A. Bare Hamiltonian 

The model we study is well-known. The system con- 
sists of ID chains, which are arranged into a square array 
to form a three-dimensional (3D) system. Its Hamilto- 
nian has the form: 

i ij 

where Hf^ is the Hamiltonian of an individual ID chain 
i. 

Hl"^ = Hf"' + HPP, (2) 
H^^- = -^VFY.P I d^--^U^^%.i)--^ (3) 

pa 

Hf = \ dx^ - J2k^'&2k^%-^-2kvi (4) 
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P2fep = Xl^R<T^La' P-2fep 

(7 



^2fep' 



c — Cit 



(5) 
(6) 

(7) 



Here index p = ±1 labels different chiralities of ID elec- 
trons, right- movers V-'Rcr (p = 1) ^^nd Icft-niovcrs ip-i,cr 
[p = —1). Vector T is composed of three Pauli matri- 
ces. The coupling constants g4.2kF £^nd J2kp are positive, 
which corresponds to repulsion between electrons. The 
model's microscopic cutoff is denoted by A. 

We expressed our Hamiltonian H^'' is a somewhat un- 
usual form. In a more traditional notation this operator 
looks as such: 



I 



dx 



52fcF = 
■^2feF = 



9lP2kFiP-2kFi + 92 ^ Pl^aiPKa'i (8) 
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Both forms are absolutely equivalent. Eq.(4) suits us 
more, for it explicitly shows the couplings of the density 
waves. 

Closely located chains are coupled by single-electron 
hopping H^""^ and electron-electron interaction H^^ : 



rrhop 

ij 



H^j = j dx\^gi{i- j)piPj 

+ 92kF {i - j) (P2fcFiP-2feFJ + H.c.) 

P = ^PPCT- 

pa 



.(11) 

(12) 

(13) 



We assume that our microscopic model is characterized 
by the following hierarchy of material constants. The 
anisotropy ratio is small: 



r = t/vFA. < 1. 



(14) 



In addition, the 2kF coupling constants are smaller than 
the coupling constant corresponding to interactions of 
smooth components of the density, and the in-chain in- 
teractions are larger (or much larger) than the inter-chain 
interactions: 



92kp ^9o < 92kF ^ ^2feF - 54 < -yp. 



(15) 



The smallness of the transverse couplings as compared 
to the in-chain coupling constants assures that at high 



energy the system may be viewed as a collection of weakly 
perturbed ID conductors. The smallness of all coupling 
constants as compared to vp indicates that weak coupling 
arguments may be applied. 



B. Effective description 

It is tempting to study the low-temperature phase di- 
agram of H, Eq. (1), with the help of the mean field ap- 
proximation. Yet, one has to keep in mind that this idea 
is wrong. It is demonstrated by Prigodin and Firsov [26] , 
who investigate the renormalization group (RG) flow of 
the QID metal, that at high energy the Cooper channel 
and the particle-hole channel are coupled, the usual lad- 
der summation is not adequate, and it is necessary to use 
the parquette approximation. Failure of the ladder ap- 
proximation implies the failure of the mean field theory, 
since the two approaches are equivalent. 

Fortunately, this coupling between the channels is a ID 
feature, which disappears at sufficiently low energy. In- 
deed, it is also proven in [26] that the weakly interacting 
QID Fermi system experience the dimensional crossover 
at low energy; below the crossover the channels decou- 
ple, and the ladder approximation (hence, the mean field 
theory) is valid again. 

Thus, if we need to know a low-energy properties of the 
model (e.g., the phase diagram), it is enough to derive 
the effective Hamiltonian valid below the dimensional 
crossover, for this Hamiltonian may be analyzed with the 
help of usual mean field approximation. 

How does this Hamiltonian look like? This question is 
addressed in several theoretical papers [1, 2, 3, 4]. These 
papers discuss in detail the dependence of the effective 
coupling constants on the bare one. However, for our 
purposes it is enough to guess general features of the 
low-energy Hamiltonian. Our conjecture is based on two 
assiimptions about RG flow: 

(i) at high energy the transverse single-electron hop- 
ping is the most relevant operator in the problem, and 
the crossover occurs when the effective transverse hop- 
ping becomes comparable to the running cutoff; 

(it) ad, high energy the SDW and CDW susceptibilities 
are the dominant; among these two the SDW suscepti- 
bility prevails. 

What do these statements mean physically? Assump- 
tion (i) is valid provided that bare interactions are suf- 
ficiently weak (for an accurate criterion of the single- 
electron hopping relevance one can consult e.g., [28]). 
It guarantees that our low-energy model is Fermi liquid 
with weak effective coupling constants and low effective 
anisotropy: 



(16) 



When (i) is violated, the system freezes into an ordered 
state (typically, SDW or CDW) before reaching the Fermi 
liquid regime. Needless to say, our analysis is inapplicable 
in such a situation. 
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Assumption (ii) is a consequence of the electron repul- 
sion combined with the fact that high-energy physics is 
purely one-dimensional. It is known that the ID metal 
has strongly divergent susceptibilities toward SDW and 
CDW. Of these two, the former is stronger due to the 
in-chain backscattcring gi . Because of all this, the SDW 
susceptibility prevails in the high-energy regime. 

This effect has nothing to do with the nesting prop- 
erties of the actual Fermi surface, which is a low-energy 
feature. Moreover, one can say that this abundance of 
the high-energy modes, enhancing SDW correlations re- 
gardless of the nesting, is a peculiarity of QID metal, 
which makes its physics so unusual. 

Taking (i) and (ii) into account we can write the fol- 
lowing effective Hamiltonian: 



The spin-spin transverse interaction term Hf-^ is equal 
to: 

= j dx[j^(i-3)Si-Si+ (18) 

■^2fep(« - i) (S2fcpi • S_2fepj + H.C.) 

S = ^w:V^.V:- (19) 

This term, although absent in the microscopic Hamilto- 
nian, appears at low energies. 



Other terms of Eq. (17), Hf"", i/f, Hf°^, and iJf/, 
have the same structure as the corresponding operators 
without tilde {Hf-^, fff, 7J^/p, and i?f/), but the former 

have renormalized constants {v-p instead of tip, t instead 
of t, g'a instead of g's). On top of this, the cutoff of 
the effective theory is much smaller than the microscopic 
cutoff: A <C A. 

Hierarchy of the effective coupling constants differs 
from Eq. (14) and Eq. (15). At the dimensional 
crossover the transverse hopping becomes comparable to 
the cutoff [sec Eq.(16)]. The effective system remains 
anisotropic (for example, its Fermi surface consists of 
two warped sheets disconnected from each other), yet, 
this anisotropy is not as strong as the anisotropy of the 
original microscopic system. 

Since at high energy the dominating fluctuations are 
SDW and CDW, 2A:f coupling constants are enhanced: 



52/cF > 54, 

.94f > aa- 



(20) 
(21) 
(22) 
(23) 



All coupling constants are smaller than the renormalized 
Fermi velocity Vy ■ 

It is tempting to declare that, since SDW correla- 
tions dominate over CDW correlations in the high-energy 



regime, the SDW coupling constant is bigger than 
the CDW constant 5^^. However, this is not necessary 
true for bare J^^^ is zero, while ^^^.^ ^ 0. This might 
affect the outcome at the crossover scale. 



III. PHASE DIAGRAM 

In this section we apply the mean-field analysis to the 
effective Hamiltonian H, Eq. (17). 

A. Density wave phases 

The low temperature phase of the effective Hamilto- 
nian depends on the nesting properties of the Fermi sur- 
face. Assume first that only nearest-neighbor hopping 
amplitude ti is non-zero. In this case the Fermi surface 
nests perfectly. The SDW susceptibility is equal to: 



1 /2€.fA 
XSDW = —r- In -— 

TTfF \ 1 



(24) 



hk^+zgikp may 



The CDW susceptibility is the same. 

As it is obvious from Eq. (18), the coupling con- 
stant for SDW is equal to gsDW = ^2feF + ^^-tk-p, 
where z is the number of the nearest neighbours. The 
usual mean-field equation for the critical temperature 
5SDwXSDw(7sDw) = 1 gives us the formula for Tsdw: 

T^^^ = 2f;FAexp[-7r5F/(J2fep + zJij,^)]. (25) 

The subscript 'max' is to remind us that at perfect nest- 
ing the transition temperature is the highest. 

The CDW coupling constant gcDW 
be larger or smaller than gsDW. depending on the bare 
values of g, g^, and J2feF- The density wave type is de- 
termined by comparison of the coupling constants: if 

5CDW = 52fep + zg2kp > 5SDW = '^2fep + zJikp, (26) 

the ground state is CDW, otherwise, it is SDW. 

When the nesting is spoiled (for example, by introduc- 
ing next- to- nearest neighbor hopping amplitude t2), the 
density wave critical temperature decreases. This hap- 
pens because antinesting destroys the divergence of the 
susceptibility. For example, one might write for SDW 
(the case of CDW is identical): 

1 f In (2iFA/T) , if T > i2, 
XSDW oc — X <^ ) ^ _ (27) 

^VF ln(25FA/i2) , if T<t2, 



(28) 



where t2 is the renormalized value of t2 ■ 

When ^2 is bigger than certain critical value: 

2 \ j-C ^max 

12 > l2 — -'SDW 



the mean-field equation 5sdwXsdw(T'sdw) = 1 has no 
solution. Thus, exponentially small t2 is enough to de- 
stroy the density wave phase. 
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B. Superconductivity 

When the antinesting destroys the density wave, the 

system becomes superconducting. To demonstrate this 
let us introduce the following set of Cooper pair creation 
operators: 



(29) 



Operator (Ay)^^' creates a Cooper pair composed of a 
left-moving electron of spin a on chain i and of a right- 
moving electron of spin a' on chain j. 

Matrix Ay may be symmetrized with respect to the 
chain indices as well: 



A«/» = - ( 
'3 2 V 



(30) 



The superscript 's' ('a') stands for 'symmetric' ('antisym- 
metric'). 

Further, it is convenient to write A*j'" as a sum of three 
symmetric matrices irr^ and one antisymmetric matrix 



IT 



y- 



d^/" • (ifr^) + A^f ir^ 



/a. 



(31) 



where r = {t^jT^, t^) is a vector composed of three Pauh 
matrices. Operator A*j^" i'i-ij"') creates a Cooper pair in 
a singlet (triplet) state. 

Using these operators we can rewrite H^j' and H^j^■. 

Y,Ht[ + Hf^ = -Y.I dx\l,._,.i:^l,{/:^l,)^ (32) 



+fl.,A^(A^)t +5/d^ . (d^)t +5^d^ • (d^)t 



+ 



where the ellipsis stand for the terms, which cannot 
be expressed as a product of a Cooper-pair creation 
and a Cooper-pair destruction operators (for example, 

V'lctV'l<t'V'l<t"'^L(7"')- The coupling constants are: 



-y2 = 

Qxy = -QJikp + 254f + ^-^0 ~ ^5^' 
~9f = 2^2fe + 252"-^ - 2 Jo^ - 2g^, 
9'f = -'^Jik^ - 2524p - 2Jo^ - 250^. 



(33) 
(34) 
(35) 
(36) 



We see from Eq. (36) that the order parameter (d") 
is always zero, since the coupling constant g'j, is always 
negative. 

Three other order parameters may be non-zero. Con- 
sider first dx2^y2-wa.ve ((A') ^ 0). This type of super- 
conductivity is at least metastable, if: 



gx2-y2 > <^ SJ^F > 92kF- 



(37) 



In the latter inequality we neglected gg and Jq- for they 
are small [see Eq. (22) and Eq. (23)]. 



Triplet /-wave superconductivity ((d*) 7^ 0) is always 
metastable, since g/ > [provided that Eq. (22) and 
Eq. (23) are satisfied]. This guarantees that, after the 
density wave is destroyed by the antinesting, the QID 
metal becomes a superconductor. 

Singlet dxy-wavc superconductivity ((A") 7^ 0) is 
metastable, if g^y > 0, which is equivalent to: 



92kF 3 . 



(38) 



The true, stable ground state is determined by com- 
parison of the mean-field transition temperatures for dif- 
ferent superconducting order parameters. These temper- 
atures are the solutions of the equations gaXa{T^) = 1, 
where a is either /, or dx2_y2, or d^y 

Let us compare first T^~y^ and T/. These two order 

parameters have identical orbital structure. Therefore, 
their susceptibilities are the same: Xx^-y'^ ~ Xf- Con- 
sequently, in order to determine the relative stability of 
dx2_y2-wave and /-wave we must compare g^^-y^ and gf. 
Specifically xf-y'' > T/, if: 



2feF 



(39) 



Eq. (39) implies that spin-density fluctuations, which en- 
hance J;^F' f'ivor dx2-y2-wa.ve over /-wave [5, 15]. Thus, 
proximity to the SDW phase promotes the former type of 
order. On the contrary, close to CDW the charge- density 
fluctuations intensify, effective coupling ^^^.^ grows, ad- 
vancing the /-wave superconductivity. 

To make this argument more concrete, consider the 
transition separating CDW and superconductivity. Sta- 
bility of CDW implies that Eq. (26) is fulflUed. This 
inequality may be rewritten as 



52/sF - J^kF > ^ {9x^-y^ - 9f) ■ 



(40) 



Observe that the bare constants satisfy [see Eq. (9) and 
Eq. (10)]: 



92kF < J2kF- 



(41) 



If we assume that the renormalized coupling constants 
g2kF and J2feF satisfy the same inequality: 



then we obtain: 



92kF < J2kF^ 



9x^-y^ - gf <0. 



(42) 



(43) 



This means that CDW cannot have common boundary 
with dx2_y2-wave superconductivity. 

However, inequality Eq. (41) is not equivalent to 
Eq. (42). Therefore, this argumentation points to a trend 
rather than establishes a hard connection between the 
density wave type and the order parameter symmetry. 

The coupling constant for d^y-wave is always smaller 
than gf. 



9f - 9xy 



&JikF > 0- 



(44) 
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This does not necessary mean that /-wave superconduc- 
tivity always overshadows the d^y-weivc superconductiv- 
ity. These two have different orbital structure (the former 
is symmetric with respect to inversion of the transverse 
coordinates, while the latter is antisymmetric). Thus, it 
is possible to choose the density of states in such a way 
that Xxy > Xf- If this happens, and if the difference 
between gf and gxy [see Eq. 44] is not too large, the dxy- 
wave order parameter may be more stable than /-wave. 
Therefore, we conclude that, when dx2-y2-wave is unsta- 
ble, both /-wave and d^y-wave can be possible choices 
for the symmetry of the order parameter. 

It is also important to note that the superconductivity 
becomes possible only in a system with sufficiently pro- 
nounced anisotropy (r <C 1). Otherwise, the ID renor- 
malization is weak, and 

92kF~92k^, 9o^9q, Jo,2k^^0- (45) 

As a result, instead of Eq. (22) and Eq. (23), one has to 
use Eq. (15). In this case all superconducting couphng 
constants are negative because they are dominated by 
—2gQ term. 

On Fig. 1 we present the phase diagram, which emerges 
from our discussion. It is drawn on the temperature- 
pressure plane. The effect of the pressure is twofold. 
First, it increases the transverse hopping, which leads 
to increase of the anisotropy ratio r. As a result, the 
dimensional crossover temperature grows, and the super- 
conducting critical temperature decreases. 

Second, the pressure increases the next-to-ncarest 
neighbor hopping, spoiling the nesting properties of the 
Fermi surface. Thus, the density wave transition temper- 
ature vanishes, when the pressure exceeds some critical 
value Pc- This explains the major features of our phase 
diagram. 

IV. DISCUSSION 

We demonstrated under rather general assumptions 
that (i) QID metal has a superconducting ground state, 
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and ( ii) the symmetry of the order parameter is sensitive 
to the details of the interaction and the density of states. 

This implies that the "universal" phase diagram of the 

Bcchgaard salts [25] is a robust feature of QID metals, 
easily reproducible theoretically. At the same time, pre- 
diction of the order parameter symmetry on the basis 
of the microscopic model is virtually impossible for it re- 
quires very accurate calculations of the competing states' 
energies. 

One must remember that different QID organic su- 
perconductors, despite many similarities they share, may 
have different symmetry of the superconducting order pa- 
rameter. Furthermore, it is possible that an individual 
sample experience a phase transition between different 
types of superconductivity, if the external parameters 
(pressure, magnetic field) change. It is possible that 
such a phenomena is indeed observed experimentally: 
within the superconducting region of (TMTSF)CIO a 
phase transition is detected when the external magnetic 
field exceeds some critical value [24]. 

Our approach, however, has some limitations. Namely, 
it is not applicable in situations, where the dimensional 
crossover occurs due to transverse interaction rather than 
transverse hopping [29] , are beyond its reach. 

To conclude we discussed the symmetry of the order 
parameter in QID metallic Fermi system. It is demon- 
strated that for a certain class of QID superconductors 
the order parameter symmetry is a non-universal feature. 
It is determined by a delicate interplay of microscopic 
constants characterizing the system. The order param- 
eter type may easily change in response to variation of 
external (pressure) or internal (doping) factors. 
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FIG. 1: Qualitative phase diagram of our model on the 
plane prossurc-tcmpcraturc. Growth of the external pressure 
acts to increase the transverse hopping, which decreases the 
anisotropy and increases the antinesting. Solid lines show 
second-order phase transitions into the density wave state 
(SDW or CDW) and the superconducting phase. Dashed line 
shows the first-order transition between the density wave and 
the superconductivity. The dotted line at high temperature 
marks location of the dimensional crossover. 



